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1 . I NTRODUCTIOTn 



In [223, Plato states that his teacher, Theodorus of Gyrene, 
was th<^ first individual to prove that VS , , . . . , Vl 7 were 
irrational 12, p. UtiJ . Since then many scholars have wondered 
how he did this with the mathematical tools he had at hand, 
and, perhaps more puzzling, why he stopped at the number 
seventeen. An answ'er w’hich is almost certainly incorrect but 
which is nonetheless intriguing was given by J. H. Anderhub in 
1918 Csee [3,5,6,1033). He suggested that Theodorus might have 
used the following construction. Form the triangle with 
vertices at CO, 03), C 1 , 03) , Cl, 13). Note that the hypotenuse has 
length V2. Starting now at the point Cl, 13), draw a line 
segment of unit length perpendicular to the hypotenuse and in 
the counterclockwise direction; this gives a new triangle with 
this line as one side, the previous hypotenuse as the other 
side, and a hypotenuse of length V3. Continuing in this way, 
a spiral figure Ccalled the q'ucLdrat-umrz^e'lschn&ck& [103, or 
square-root-snail J emerges with the numbers Vn as the lengths 
of the radii CFig. 1 3> . The last hypotenuse drawn before the 
figure overlaps itself has length Vl7, and this, it was 
suggested, was why Theodorus stopped at seventeen Csee C 53 3). 

Although this is almost certainly the wrong answer, the 
sequence of points formed by this process was studied on its 
own merits by Hlawka CE103; see also [2433). He proved a 
number cf results concerning this construction. The study was 
then pirKed up by Davis [33, who noted that the points could 



Cl . ID 



i xii ^ y\. ^xcxiiS' ^ ^ x«w'C7x c^UXs^lx 

Z = Z + vZ / I Z I 
r»-*-l n n r> 

Cwi'th 7:^=1 D and that, they we:- e asymptotically an Archimedean 
spiral. He then interpolated an analytic curve to these 
points, which he called the Spiral of Theodor us. 

Noting that Cl. ID could be considered Euler’s method for the 
ordinary differential equation 

z = z/ I z I 

with unit step size, Davis suggested as a generalization of 
Cl. ID the iteration 



z = oz + /?z / I z I 

n-*-i n n n 

with a,/?eC and a given complex number. This iteration, 

called the Complex Generalized Theodor us Iteration CCGTID, is 
studied in [3,15,203. Davis later suggested the further 

gener al i zati on 

V = A^V + B^V xlIV II Cl. 2D 

n-^l n r> n 

where A and B are real mxm matrices, is a given nonzero 

m-vector , and 11*11 is the Euclidean vector norm. This 
iteration, named the Generalized Theodorus Iteration COilD by 
the present author, displays a great variety of strange 
attractors for appropriate choices of A and B CFigs. 2-5D , and 
includes the complex GTI as a special case; it is studied in 
[4,15,13,193. In this report will summarize some of the 
inf orma <,.ion contained in [163 and [193 Cadding some material 
from [153D with respect to sufficient boundedness conditions 
and asymptotic estimates for this iteration. Since we are 
interested in boundedness and asymptotics, in stating our 



result? we will ignore the case where V is such that V is 

O n 

the null vector for some n Cin which case the iteration 
stops!) . 



Boundedness 



Let us try to find sufficient conditions for boundedness of 
the GTT . Taking norms on both sides of C 1 . 2!) gives 



II V II < HAH • II V II + HBIl 

n-t-l n 

Iterating this relation, we find that 

n- 1 

IIV II < llAll'^-IIV II + IIBII-r hah' C2.1D 

n O 

t =o 

and clearly if II A IK 1 then in the limit the norm is bounded by 

IIBII/Cl - \\A\\:> 

Hence lAIKl is sufficient for boundedness of the orbit C f or 
any • From this and the similarity of Cl . 2!) to the linear 



system 



V = A^V CP..P:0 

n 

which has a bounded solution for every when the spectral 

radius of A is less than unity Cand at least one unbounded 
solution when the spectral radius of A is greater than uni tyj , 
it seems reasonable to conjecture that pC AD< 1 suffices for the 
boundedness of C 1 . 2!) , where pC AD is the spectral radius of A. 
That this is in fact the case follows from two lemmas. We 
omit th‘? rather technical proof of Lemma 1. It is not hard to 
establish this result by using an approach similar to that 
LHwCi in fiV0 SLnci six but it is rather tedious to 

carry cut the proof in detail Csee [153D. Lemma £ is a 
discrete version of an ordinary differential equations result 



.•1 



Lemma 1 • 



It I! V II is unbounded Lhen IIV II — >- 00 . 

n n 

Lemma S: Lei A be a mairix wi ih pC A!)< 1 and suppose ihat BC n!) 

is a sequence of matrices with lIBCnDll — as n — >oo. Then all 

sol uti ons of 

2 =CA + BCn:):)‘2 

n-»-l n 

tend to zero as n — ►ao. 



Proof: Since pCA3<l there exists a matrix norm I! • II such that 

A 

IIAll <1. We will choose this to be a natural matrix norm and 

A 



will 


al so denote 


by II • II the 

A 


vector 


norm 


whi ch 


i nduces 


i t . 


Let 


e equal 


Cl- 


II A 1! D. NoVf’ 

A 


si nee 


all 


ma t r i x 


nor ms 


are 


equi 


val ent i n 


t he 


sense that if 


II MC n'J II 


—►0 

Ot 


t hen II MC 1! ^ — ►O , 


we 



have from ovir hypothesis that lIBCnDll^ — >0, Hence we can choose 



an W such that^ II II 1 c^r all iilNt 

A 



Now consider iterating 



- C A BCnJJ '7 

n-H n 



C f or some z D N times to produce the vector z Set w =z and 

O NOW 



= C A + BC N + nD ‘ O' 

n**-! n 

Vn^O. CJ early to =z for all i>:0. Now 

\ N+i 



o = C A + BC N + nD J X • ■ • A + BC Nj D xto 

T>^± O 



and so 



llo) I! < IIA + BCN H- nDll ^ — ^IIA h- BCNDII ^IIu' II 

n-fl A A A O A 

n 

= llw II n llA + BCN + iDII 

O A * * A 

i = 0 



Tnen 

n 

llo) II < llco II n CIIAII + liBCN + iDll :> 

nH A O A ** A A 

V =0 



< llco II n - eD + e/S:) 

O A * * 

V =0 

r> 

= llu' II n Q1 - 

O A • * 






which clearly tends to zero as n — >oo. Hence llco I! — ►O as n — 

n A 

SO that Hz li — >0 as n — >oo, and so by the equivalence of vector 

n A 

norms on [R^ we have that liz II — ►O as n — ►oo. This completes the 

pr oof . □ 

We may now show the desired result. 

Theorem 1: If pQA'^<l then the solution of C 1 . 2D is bounded. 

Proof: Suppose the contrary and choose a V such that IIV II is 

O n 

unbounded. By Lemma 1, IIV II — ► 05 . Define BCnD=B/'IIV II. By the 

r» r» 

hypothesis, lIBCnDll — >0. Consider the iteration 

z = CA + BCnDD^z 

n-4-l n 

with z =V . Clearly z =V for all n. The conditions of Lemma 
c» o r. o 

2 are met and so IIV II — ►O as n — ► od. But this contradicts IIV II 

n n 

unbounded, which establishes the result. □ 

It i S: also possible to show that if pCAD>l then there exists 
a V such thiax. IIV II is unbounded Csee section 5D . In this 

O n 

sense the GTi is sur pjr 1 si ngl y similar to the linear system 

C 2 . 2D . 

In tt“ particular case A=0 a detailed analysis is possible 
Cl6,17] 7>ji s case is related to t.he power method for the 

numerical soltjtion of eigenvalue problems. 

3. Asymptotics 

This leaves the question of boundedness for the case pC AD =1 . 
It IS possible to construct examples showing that when A is 
uni spec T.ral the orbits of the GTI may be always bounded for 
one choice of B and always unbounded for another choice of B. 
However, a general statement about the asymptotics of the 



R 



i ib^i A 



4 4 £=:'♦-• r 4 -j WD-»^ HiX' A i >r ^ 4 " Kf^it H liitH 4 i ^ 

J. tAS^ l S.^11 **K_4. ^ illCAS.^ 'v- f-i jL^-> Oa SUiiS^VwO ^iiCAV 

pCA!)<l and any eigenvalues of A ‘that, lie on t.he unit circle 
are simple Ci.e. have as many linearly independent 
eigenvectors as their mul ti pi i ci tyJ . In this case the linear 
system CZ.Z'J has the property that all solutions are bounded 
Cal though in general the bound will depend on • Before 

stating the corresponding resul-t. for 4. he GTI , we state: 
Hukuwara's Theorem: Let z be any solution of 



z = CA 4 Ect:):)^z 

t+i t 

for t=0, where A is a constant matrix and 

CC' 



C3. ID 



J] llBCtDll < cn 
t = 1 

Suppose that. A is of bounded type, so that all solutions of 
y =A^y are bounded for t>0. Then every solution of C3.1D 
is bounded. 

TTii s result, an analogue of a theorem from differential 
equations, may be found in [£13. The proof there uses the 
vector norm, but the argument is valid for any vector norm and 
compatible matrix norm. We will now use it to prove an 
asymptrtir resuit for the GTI. 

Theorem £: Let A be of bounded type in C1.£D. TTien 

CX) 

r i/'iiv II =00 

^ n 

r> = O 



Proof: Suppose the contrary, i . e. that for some we have 

CD 

T; 1 /II V II < 00 C3. 21) 

n 

rt = O 



and construct such a sequence <V > by using this V and Cl . 2D . 

r> O 

IDefine BCnD=B/llV II and consider the iteration 

rt 

Y = CA + BCnDD^y 

n*l n 

where set Y =V . Clearly Y =V for all niO; the two 

O O ri ri 



Now 



sequencers are identical. 

00 00 
T. llBCn:)ll = 5] IIBII/IIV II 

n 

n = O n = O 

00 

= iiBii r i/iiv II 

n 

n = O 

which by supposition is bounded. Therefore by Hukuwara’s 



Theoren., IIY II is bounded. Then IIV II is bounded, that is, 3M>0 



such that IIV IK M for all n. But then lll/V ll>l/M and since 

n n 

i/IIV II— x->0 we cannot have C3.BD. This contradiction 

n 

establishes the result. □ 



III fact a slightly more general result holds. If 11*11 is any 



vector norm and O is any function that is continuous on 



rO,ao!) and positive on C0,oo!), and if A is of bounded type, then 



any solut-ion of 



V = A^V + x^CllV ID 

n-*>l n n n 



must be such t-hat 

00 

2 1 IIV ID =00 

n = O 



Thii s ma ^ be shown in a completely anal ogovis manner Csee [153 



for details and applications!). 



Thieorem £ establishes that, w'hen pCAJ=l but A is of bounded 



type, if the orbits diverge then they do so in a way that is 



not. much worse than linear Csince the sum of the reciprocals 



of the noT“TTi'= diverges^). In fact if IIAII=1 in the spectral norm 



then it is a fact that the divergence is at most linear, as 



can be seen by setting IIAII=1 in C2.1!) to get 



IIV II < IIV II + nllBlI 

O 



as the f symptoti c estimate. We will later consider this point 



in grea*..er detail when we show in section 5 that if A is of 
bounded type then the divergence is always at most linear. 



4. Green's Function Representation 



In order to refine our asymptotic estimates of the previous 



section, we will need a formula for the solution of the 



iteration Cl . 2D . Given an iteration of the form 



z = ACtDxz + o) 

t t-i t-i 



with 0 .' and z given vectors and ACtD a time-varying matrix, 
to ^ ^ 

the solut.ion can always be w'ritten in the form 

t 

z = r Y -Y'^ w + Y -Y"*-c 

1 + 1 is S t + 1 o 

s = o 



where c is the initial condition z , and <Y > is a fundamental 

o t 

matrix set of solutions for 



z = AC tD^z C4. ID 

1 t-i 

Cthat is, Y =I and Y =ACnD^Y D [211. For the GTT , the 

O n+l n 

matri X A is constant, and z =V is given. Furthermore, Y =A^ 

o o ^ n 

for the linear system C2.2D corresponding to C4.1J, and 



(o = B^V /IIV II 

i t t 

so that 



V = + r /iiv II C4.a:> 

n +1 o ^ t i 

1=0 

is the solution of Cl. ID. Thiis is called the Green’s Function 



Representation of the solution. We will use this form of the 



GTI to address the questions of boundedness and asymptotics. 



5. A Direct Approach to Boundedness 



The approach used in the second section to show that pCAD<l 



suffices for the boundedness of Cl. 2D is the one originally 



G>rupl(z>y^d by ih& author in [ 1 5J and also appears in fl8]. A 
more direct method, however, was used by the author in [193, 
and we now give this method. 

Suppose that pCAD<l. Then there exists some natural matrix 

norm II • U such that II All < 1 . Using II • II to represent the 
o a a ^ 

vector norm which induces thi mAi.rix norm as well, we have 

from C4.23) that 

n 

r /IIV II 

^ i t II a 

= o 

^n>0. Thus 

n 

II V II < II AlO*‘ • II V II + r IIAII* • IIBII -k 
n +1 a o o a ^ a o 

t = o 

where t is a positive constant such that 

II VII /II VII < k 
a 

for all nonzero vectors V. Such a k exists since all vector 
norms are equivalent on [123. Summing the geometric series 

and col lee- ting like terms gives 

IIV II < eCnD + k-llBlI /Cl - llAll j 

o< a OI 

w'here cv. ri.; >o as n — >a-). Hence IIV II is asymptotically bounded, 

n Ot ^ r- 

and by the equivalence of vector norms so is IIV II. An 

r> 

asymptotic bound for IIV II is given by 

n O* 

k'llBII /Cl - II All J 

a O' 

in terms of the o-norm. When the spectral norm of A is less 
than unity the corresponding bound using that norm, 

IIBII/Cl - IIAIID, 
is generally fairly tight. 

We have shown that pCAJ<l suffices for boundedness. In an 
anal ago. js way we can show that if pCAI)>l then there must exist 
some V such that IIV II is unbounded. CA similar approach has 

o n 



IIV II < 

rt+1 Ok 



ia""‘V 



10 





. >- >-4 


U-,* V 
V. 


. 4 . >-,i IsJCAl 1 ^ 


w. r 'I e, 1 ■'1 

lj-Oj-.*’. rOi 


> 


cronsi deri ng 


C4. 2D , 


we have 


that 














IIV 

n-n 


II > 


II a''** V II 

O 


- 1 


1 E a''"Vb«v^/iiv II j 
1=0 






C usi ng 


■che spectral 


nor mD . 


Let X be an eigenvalue of A such 


that 1 


X j =pC AD 


and 


choose 


V 

o 


to be an eigenvector of A 


associated with 


the 


eigenvalue X. Then we have 








IIV 

n+1 


II > 


iix'^^V 11 

o 


- 


r> 

E /IIV II 

i 1 

1 =o 







and since the norm of the sum is at most the sum of the nor ms > 

n 



IIV II > 



IX C* IIV 



- E iiAir • IIBII 



1=0 



p IXr""* - IIBII -Cl - llAli'^^b/Cl - IIAO 

O 



wi t h p = li V' II 
o o 



then c>0 and 



I f we 1 et 

c = IIBII/'C II All - ID 



n+l 



p |X - cllAlC'^ + c 

o 



in term? of c. Now pC AD = | X | -t-e for some e>0, so 



liV II > 

n-^l 



1 ^ 

Pol^l 


-cC|X| +eD -^'C 








ixr"‘cp 


- cCl + €/|X +c 

o 


C5, 


. ID 


we may 


now choose a to be 


gr eater 


than 



c in order to have IIV II — ><x>. Otherwise, the term involving e 

will grow' arbitrarily large, again giving II V II — >oo. 

r» 

Again, w'e remark that the GTI Cl. ID has boundedness 
properties similar to those of the linear system C2.2D. In 
the next section we will derive improved asymptotics for the 
indeterminate case where A has unit spectral radius. Before 
doing so let us collect the information above as a theorem. 
Theorem 3: Consider the iteration Cl. 2D. If pCAD<l then the 



iteraticn is bounded for any V^, and if pCAD>l then there 



11 



&XJ. sis 



udi cJ^at the iteration is unbounded. 



a V 

o 



sTi. Improved Asymptotic Estimates 



From C5.1J it is immediately apparent that when pCA!)>l the 
divergence will, in general, be exponential. Tbi s again 

leaves the case pCAD=l. As before, let us first look at the 
case where A is of bounded type. Using the spectral matrix 
norm, C4.2J gives 

n 

C6. l!) 



II V li < II + 

n-H O 



r ‘xBxV /II V I 

t t 

t =o 



< K r HA 

1 

t =o 

where K >0 is a bound on 

1 



IBI 






C which in general will 



depend on • Now since A is of bounded type, II A 

boundeo , that is, IIA^II<K for some K>0. Thus 



is also 



IV II < K + IIBli T K 

n-M i ^ 

1=0 



< K + K n 
1 2 

showing that the divergence is no worse than linear. 

If pCAJ-1 but. A is not of bounded type then it is possible 
to show that II A^ II =<DC n^3) Cin fact II a’^ II =€>C n^ where p<m is the 
maximal degree of all nonlinear divisors of A associated with 
uni modular eigenvalues!), where A is mxm. Hence from C5. 1!) we 
have 

n 

iiv II < iia’^**ii • iiv II + r iia‘ii-iibii 

n+l O 

t =o 

and clearly IIV II diverges at most as . CIn fact, it 

diverges as €>C n^!) , and hence in the worst case as ^Cn^!)!). 
Hence, when pC A"J =1 , if the iteration di varg^s then it does so 



12 



i ^ •'^1 


J 1 


•'f i 'i 1 1 

C'. .j. a. j> ^ 


Wc::- -ilci 4Hi>r 's»r 

• IS- J J. S_- S.< i i X ^ O*. O'. 


4 

•w 1 Jill 


Theorem 


4: 


Consi der 


Cl . S j and suppose 


that pC Aj =1 and the 



iteration is unbounde-d. If all ©i g<anval u»a3 of A wi bh unit 
modulus^ are simple then the divergence is no worse than 
linear. Otherwise the divergence is no worse than polynomial 
of order m, where mis the order of A. 

7. Li mi t Sets 

Let us define FCA,B,V i to be the limit set Cor a.^-1 i mi t. sei. 

o 

L13]J of the iteration; that is, YerCA,B,V D if there exists 

o 

an i ncT-f==.'?ri r..n subsequence nC i such that V — ►Y as i — ►od when 

n<L> 

the in.itiaJ condition is V . Two theorems follow easily from 

o 

the material in [133, 

Theorem 5: Suppose that 

bounded. TTien rCA.B,V 3) 

o 

V — ^rCA.B,V j as n— 

n O 

Proof: ClearJy the closure of F is simply F and so F is 

closed. The boundedness of <V > implies that F is also 

n 

bounded. Hence F is compact. C by def i ni ti on3> . By the 

Bol zanc>-Wei er str ass Theorem [S33, F is not empty. 

Now V — >r means that infCllV -YlliYeFJ — as n — ►co Cl 33. 

ri n 

Suppose the contrary. Then we can find a subsequence of <V >, 
indexed by, say, nCiD, which remains a finite distance from F. 
Now V is a bounded infinite sequence and so by the 

Bol zano~Wei er str ass Theorem it contains a limit point. But 
that po: nt must be in F, which is a contradiction, d 



lor s om«rr t h «=r 

o 

is a nori -empt y 



sequtrric ^ 


j' \r \ 

V. » -* 

n 


i s 


compact 


set 


and 
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Theorem 6; 



Suppose that for some the sequence <V > is 

boundt=' ■ ar.d 0^rCA,B,V . Then rCA,B,V J is positively 

i nvar i ant. . 

Proof: E>ef i ne 

TCV!) = A^V B^VXIIVII C7.1D 



so that 



V = TCV D 

n-»*l n 

is the iteration under consi der ati on . Since Oi^T, the 

tr ansf cr ma ti on T is continuous in some neighborhood of any 
point Yer. Fix Y and choose a subsequence nCiD such that 
V — >Y . By the continuity of T, we have that 

n<i > 

TCV :>-^tcy:> 

n(v> 



as i — >cr Bui 



so that 

wher e mC i =nC i D +1 . 
invariant, n 



TCV :) = V 

r*<v> 



V —*TC YD 

rr><v> 

Hence TCYDeF, and P is positively 



In fa t, .it js clear from the above that TCYj^F whenever YeP 
and TCYl/ is defined. In other words, P\<0> is positively 
invariant. Further topological results on discrete dynamical 
systems that can be extended to this case Cwith some care 
taken near the origin!) may be found in C143. 

Suppose that YePCA,B,V D, so that V — ►Y for some 

O n<v> 

increasing subsequence nC i !) . Since the transformation 

T: [R^ — of C7,lj has the property that 

TC-\0 = “TCVj 
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it follows that the sequence generated by is -V which 

tends to ~Y ; that is, -yerCA,B,-V^ 3. In a similar way, it is 

o 

easy to s ee t hct t i A , “*0 , J and i l A , B , ^ . 

o o 

In the case of a global attractor, rCA,B,V D is independent 

o 

of V for the particular A and B under consideration. In this 
o 

case, it follow's from the above that the attractor is 
centrally symmetric, i . e. YeP implies -YeP, and that the same 
attractor is found w'hen A and B are replaced by -A and -B, 
respecti vel y. Some of the attractors displayed by this 

iteration for various choices of A and B are shown in the 
ligurtr"^', Stir'fcz:' clIso* 15j. Iserles has found an anal yt i cal 

expression for the attractors for a certain set of cases of 
the form A=oiB, cx>0 Csee C4]J. For the case A=0 the form of 
the attractor is also known [15,16,17]. 



S. Discussion 



The Generalized Theodor us Iteration is in some W'ays very 
similar to the simple linear system C Theorem 3D and yet it 
contains some very complicated structures CFigs. 2-5D . It is 
likely that the attractors are, in general, strange but not 
al W'ays chaotic Csee E9]D; some com^putati onal evidence to this 
effect has been provided by James Heyman of the Naval 
Postgraduate School in his study of the use of chaotic 
discrete dynamical systems, particularly the H4-non attractor, 
ciS pScr idorandom numiber generators C personal corumuni cati onD . 



^ cr. 
± 



, I 4 .-., — • r X w..- m~r ; x ui ^ x ^r-x 

K CJi i <a]-^^-i- J. ‘^ca J. V^i <^X S^iiCT X J. -L i i J. V^CLlJV^ 

are in progress. Due to the present plethora of 

definitions of the terms **chaos** and **f ractal '* C see [15]D it 
is difficult to prove anything substantati ve about the chaotic 
nature of the attractors but, as noted in Cll], many 
wel 1 kr nown s t r ange a 1 1 r act or s are still a wai t i ng anal yt i c al 
verification of the observed qualitative features of the 
underlying iterations and differential systems. One aspect of 
the GTJ which we feel is particularly relevant to the field is 
the vast array of geometrically distinct attractors to be 
found in the" GTI for different choices of the matrices A and 
B; even though these two matrices essentially represent eight 
parameters. the number of distinct attractors is still 
sur pr i SI ng > and we ar e not awar e of a.nother di scr ete pi anar 
map wdth this property. Further investigation of this 

pecul i a r i i 2y" seems warranted. Hopefully the figures cited here 
and in [151 will motivate research on the question of the 
number of distinct attractors which can belong to a single 
iteration, as the ’works by Gleick [7], E>evaney [S3, and 
Wiggins [S5,S6 3 have served to help motivate and inf or mi this 
author . 
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